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In previous work on adsorbate-induced surface core level shifts (SCLSs), the effects caused by 
O atom adsorption on Rh(lll) and Ru(OOOl) were found to be additive: the measured shifts for 
first layer Ru atoms depended linearly on the number of directly coordinated O atoms. Density- 
functional theory calculations quantitatively reproduced this effect, allowed separation of initial and 
final state contributions, and provided an explanation in terms of a roughly constant charge transfer 
per O atom. We have now conducted similar measurements and calculations for three well-defined 
adsorbate and coadsorbate layers containing O and H atoms: (1 x 1)-H, (2 x 2)-(0-|-H), and (2 x 2)- 
(0-I-3H) on Ru(OOOl). As H is stabilized in fee sites in the prior two structures and in hep sites in the 
latter, this enables us to not only study coverage and coadsorption effects on the adsorbate-induced 
SCLSs, but also the sensitivity to similar adsorption sites. Remarkably good agreement is obtained 
between experiment and calculations for the energies and geometries of the layers, as well as for all 
aspects of the SCLS values. The additivity of the next-neighbor adsorbate-induced SCLSs is found 
to prevail even for the coadsorbate structures. While this confirms the suggested use of SCLSs as 
fingerprints of the adsorbate configuration, their sensitivity is further demonstrated by the slightly 
different shifts unambiguously determined for H adsorption in either fee or hep hollow sites. 



I. INTRODUCTION 

Core electron binding energies of solids are sensitive 
probes of the local atomic and geometric environment 
and as such markedly changed at surfaces. A quantita- 
tive understanding of the corresponding shifts, in partic- 
ular of surface atoms coordinated to atomic or molecu- 
lar adsorbates, provides a wealth of information e.g. on 
the nature of the chemical bond [1| and may even be 
used for the interpretation and prediction of surface re- 
actions [5] . Experimentally, the shifts can be accessed by 
high- resolution X-ray photoelectron spectroscopy (XPS), 
as part of the so-called surface core level shifts (SCLSs) 
Apart from the actual so-called initial-state shifts 
due to the changed electron density viz. electrostatic 
potential at the location of the core electrons, the mea- 
sured SCLSs also contain so-called final-state contribu- 
tions from variations in the screening of the core hole 
created in the XPS process, which is influenced by the 
surface or nearby adsorbates as well. At present, a parti- 
tioning of the experimentally only accessible total shifts 
into initial- and final-state effects is only possible by 
first-principles electronic structure calculations. In par- 
ticular for metal surfaces, approaches based on density- 
functional theory (DFT) have recently proven to be very 
successful. For instance, remarkably quantitative agree- 
ment could be achieved with the SCLS measurements 
of the first-layer atoms of two close-packed transition 



metal surfaces, Rh(lll) [I and Ru(OOOl) and of their 
modification by varying amounts of adsorbed O atoms 
in well-defined structures. Moreover, these calculations 
were then able to differentiate charge density and screen- 
ing contributions and to analyse the respective physical 
basis in terms of the structures and the electronic prop- 
erties of the systems studied. 

Especially for adsorbate systems the high-resolution 
XPS spectra generally exhibit multi-peaked structures 
arising from different environments of the surface sub- 
strate atoms (first, second, or deeper layers; neighbors 
to adsorbates of different type and number). Some as- 
signment can often be reached by suitable spectroscopic 
means or by systematic variation and comparison of var- 
ious sample conditions e.g. in form of different adsorbate 
coverages. However, a comprehensive interpretation re- 
quires detailed DFT calculations which in turn necessar- 
ily need to be based on a microscopic structural model. 
In view of a potential application to less well character- 
ized systems more general relations enabling a direct data 
interpretation are therefore of interest. In this respect, 
the mentioned SCLS work at Rh(lll) 3 and Ru(OOOl) 
[1] identified an intriguing scaling in that the 0-induced 
SCLSs of the first-layer metal atoms turned out to be 
linearly proportional to the number of directly coordi- 
nated O atoms to good accuracy. With the induced SCLS 
containing only a small and not much coverage depen- 
dent final-state contribution, the source of this additiv- 



ity could be traced back to a roughly constant amount 
of charge that is withdrawn from the surface atoms by 
each directly coordinated O neighbor With ensu- 
ing work on Rh(lOO) confirming this then termed "ad- 
ditivity rule" , it was suggested that the latter may be 
a universal property of SCLSs which could be used for 
the fingerprinting of adsorbate populations in catalytic 
applications 

The cases mentioned all concern O atoms, a species 
with high electronegativity. SCLS results have also been 
obtained for H atoms on Rh(lll), (100), and (110), in- 
cluding coverage dependencies [5] , but no results exist for 
mixed layers. It is therefore very interesting to examine 
whether such an additivity rule applies also to combina- 
tions of different adatoms. To this end we have selected H 
atom adsorption and O + H coadsorption on Ru(0001). 
This choice has the advantage that not only O 0, 0] and 
H atomic overlayers, but also two well-defined (0-l-H) 
coadsorbate systems [1, |^ have already been character- 
ized in detail in terms of adsorbate phases and geome- 
tries. On this basis we have measured the SCLS spectra 
for these systems and combined them with DFT calcu- 
lations for ground state energies, geometries, electronic 
structure, and SCLS values for the various distinguish- 
able Ru surface atoms. We find excellent agreement be- 
tween experiment and theory both for the geometries of 
the stable structures and for the SCLSs of all surface 
atoms. The central result we obtain is that we find the 
additivity rule to apply quite well also for H adsorbate 
atoms, and even for coadsorbate structures containing 
the strongly charge-transfer inducing O atom and the 
roughly neutral H atom. 

The paper is structured in the following way. An ex- 
perimental section briefly summarizes the technique, the 
studied layers and their preparation, and displays the ac- 
cumulated XPS data and the derived shift values. The 
fitting procedures to extract the latter from the mea- 
sured spectra are briefly described in an Appendix. The 
next section details the computational procedures and 
compares the obtained stable geometries in detail with 
existing experimental data. We then compare the total 
calculated shifts to the experimental ones, and decom- 
pose these total SCLSs into initial- and final-state con- 
tributions. Concluding with a thorough discussion and 
analysis of the physical properties behind the observed 
"additivity rule", we can trace the latter again back to 
an almost constant — albeit compared to O neighbors 
much smaller — amount of charge transferred to each 
directly coordinated H neighbor. 



II. EXPERIMENT 

A. Structures: Geometry and preparation 

The experimental procedures were similar to those of 
ref. The SuperESCA beamline of the synchrotron 

radiation source ELETTRA in Trieste was used which 
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(2x2)-0 (2x2)-(0+3H) 




FIG. 1: Schematic top views of the examined surface struc- 
tures prepared starting from either the clean or (2 x 2)-0 cov- 
ered Ru(OOOl) surface. The color-coding shows the different 
next-neighbor surroundings of the surface Ru atoms. 



works at base pressures in the low 10^^*^ mbar region and 
contains preparation facilities and an electron spectrom- 
eter with multichannel detection. The RufOOOl) surface 
was prepared by estabhshed procedures [1,0,0,Q; clean- 
liness and order of the surface were checked by XPS 
and low-energy electron diffraction (LEED). The follow- 
ing adsorbate structures have been investigated: (1 x 1)- 
H, (2 X 2)-0, (2 X 2)-(0-hH), and (2 x 2)-(0-|-3H). For 
all structures, the adsorption sites and detailed geome- 
tries are known frornprevious quantitative IV-LEED or 
VLEED analysis [1, U . Figure [T] compiles the structures 
and the sequence of preparation. In (1 x 1)-H all hydro- 
gen atoms sit in fee sites. To prepare this saturated layer, 
the well-cleaned surface has to be exposed to high doses 
(several 1000 Langmuir) of molecular hydrogen, because 
the last few percent to saturation coverage adsorb very 
slowly jTo| . Low background pressure when hydrogen 
is introduced is important: otherwise the layer becomes 
contaminated with water [1Q|. For the two 0-l-H coad- 
sorbate structures, the starting point is a well ordered 
(2 X 2)-0 structure with O adsorbed in hep sites. The 
latter can be produced by slight overdosing with oxygen 
at 300 K to produce an O coverage of about 0.3 mono- 
layer, heating to 1250 K, then heating to 400 K in 10"^ 
mbar hydrogen for about 10 minutes, and finally heating 
to 400 K for ordering [ll[ . Saturating this structure with 
hydrogen at 90 K, followed by short anneal at 100-110 
K for ordering, produces the (2 x 2)-(0+3H) structure. 
Upon heating the latter to 155 K with 2 K/s, or faster 
to 160 K, two thirds of the hydrogen is desorbed and 
the (2 x 2)-(0+H) structure remains, as was proven by 
thermal desorption spectra (TPD) [8| and quantitative 
LEED [9] . Quantitative LEED has _also corroborated ex- 
pectations from quantitative TPD [S] that in the (0-I-3H) 
structure all H atoms sit in hep sites, while in the (O-fH) 
all H atoms sit in fee sites [9] . This site selectivity in the 
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FIG. 2: High-resolution XPS spectra of the Ru 8^5/2 core 
level of clean Ru(OOOl), measured (from top) at 333, 345, and 
400 eV photon energy. The main components are marked 
(Si: first layer Ru atoms; S2: second layer Ru atoms; bulk 
Ru atoms). The bulk value is used as zero reference for the 
surface core level shifts (top a;-axis). 



two structures has been ascribed to repulsive 0-H lateral 
interactions coupled with the tendency to maximize the 
coverage and we come back to this point below. 

The interesting gist for the SCLS measurements is that 
the various structures contain Ru first layer atoms in a 
variety of surroundings as also schematized in Fig. [T] 
with 3 Hfcc neighbors in the (1 x 1)-H layer; with 1 O 
or without adsorbate neighbor in the (2 x 2)-0; with 1 O 
+ 2 Hhcp or 3 Hhcp in the (2 x 2)-(0+3H); and with 1 O 
+ 1 Hfcc or without adsorbate neighbor in the (2 x 2)- 
(0+H). Together with the structures containing zero- 
to threefold 0-coordinated Ru surface atoms studied in 
ref. [1], this variety offers a good basis to systematically 
check on the SCLSs additivity with and without coadsor- 
bates, as well as with H adsorbates in slightly different 
hollow sites. 



B. SCLS measurements and analysis 

The SCLS experiments were performed at three dif- 
ferent photon energies (333, 345, and 400 eV) at a pho- 
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FIG. 3: Ru 3(^5/2 SCLS spectra for all studied structures. 
Additionally shown is the decomposition into the different 
peaks, where the color-coding indicates the assignment to the 
different surface atoms as explained in Fig. [1] 



ton beam incidence angle of 70°, which in the present 
setup leads to electron emission into the surface normal. 
Through spectroscopic means and detailed DFT calcula- 
tions ref. y has already estabhshed the assignment of 
the peaks to the different Ru atoms (first layer, second 
layer, bulk; modification by O neighbors) for the clean 
and the (2 x 2)-0 surface, thereby also facilitating the de- 
tection and analysis of new peaks in the H (co)adsorbate 
structures. As illustrated by Fig. [5] for clean Ru(GOGl) 
the data obtained at 333 eV photon energy appears most 
surface-sensitive in that both the first layer (Si) and sec- 
ond layer (S2) surface peaks dominate over the bulk in- 
tensity. All SCLS data and spectra shown in the follow- 
ing therefore correspond to this photon energy. The data 
acquired for the other energies has been analyzed as well, 
without obtaining significantly different findings to those 
detailed below. 

Figure[3]collects the SCLS measurements of all surfaces 
studied and shows the components derived from the fits; 
the numerical values are given in Table 3 (vide infra). 
The fitting procedures used to distinguish and disentan- 
gle the different components in the spectra are identical 
to those discussed in ref. with further details given 
in the Appendix. The best-fit parameters (asymmetry 
a, Lorentzian width L and Gaussian width G) deter- 
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mined for the clean and (2 x 2)-0 surface are close to 
those published in ref. small changes are explainable 
by the higher surface intensities obtained in the present 
measurements. In fact, we believe that the present values 
for the surface peaks are the more reliable ones. In any 
case, the central point is that all re-measured SCLSs for 
the clean and the (2 x 2)-0 surface are within error bars 
the same as in the previous work. The best-fit parame- 
ters determined for the H-containing structures are again 
similar, with further details given in the Appendix. One 
notable difference is that the components corresponding 
to Ru atoms with H neighbors had to be given some- 
what larger widths, which we interpret by inhomogeneous 
broadening due to a less than perfect ordering achieved 
for these structures. 



III. THEORY 
A. Computational setup and SCLS calculations 



was employed. 

The chosen computational approach is thus highly sim- 
ilar, if not superior to the one used already in ref. f^, 
with the notable exception that the now computationally 
feasible thicker slabs should enable a better description 
of the SCLSs of second layer Ru atoms. Recomputing 
all SCLSs of all oxygen-containing structures reported 
in ref. Q with the new computational setup, it is thus 
not surprising that we obtain with one exception val- 
ues that are identical to the published ones to within 30 
meV, thereby simply reconfirming the stringent conver- 
gence tests performed for the previous work. As expected 
the sole disagreement concerns a second layer shift in the 
(2 X l)-0 structure, denoted as S2(0) in ref. [ij, that was 
previously computed as —21 meV and is now obtained as 
+61meV. This does not affect the published comparison 
to the experimental data though, as the latter peak was 
not resolved in the measurements. 



The DFT calculations were performed using the 
full-potential augmented plane wave plus local orbital 
(L)APW-flo method [H [H] as implemented in the 
WIEN2k code The exchange-correlation (xc) en- 

ergy is treated with the generalized gradient approxima- 
tion (GGA) using the PBE xc-functional The pro- 
cedure to obtain initial- and final-state contributions to 
the SCLSs is exactly the same as detailed in ref. in 
which the initial-state shifts are simply given by Kohn- 
Sham eigenvalue differences and the final-state screening 
is computed within the Slater- Janak transition state ap- 
proach by an impurity calculation with half an electron 
removed from the core state considered. 

All surfaces are modeled with supercell geometries, 
using inversion symmetric slabs consisting of eight 
Ru(OOOl) layers with O or H adsorption on both sides. 
The vacuum between consecutive slabs is at least 11 A. 
Within the employed (2 x 2) surface unit-cells the po- 
sitions of all adsorbates and Ru atoms in the topmost 
three surface layers were fully relaxed for all structures 
considered. The gas-phase calculations for atomic H and 
O required for the binding energies are done in rectangu- 
lar supercells with side lengths (13 x 14 x 15)bohr. The 
obtained total energies are then corrected by the gas- 
phase O2 and II2 binding energies to obtain adsor- 
bate binding energies with respect to the molecular refer- 
ences. The basis set parameters are briefly summarized 
as follows: The muffin-tin spheres for Ru, O and H are 
= 2.0bohr, R^^. = l.Obohr, and R^j, = 0.8bohr, 
respectively. Inside the muffin-tins the wave functions 
are expanded up to l'^^^ = 12 and the potential up to 
^max = 6. The energy cutoff for the expansion of the 
wave function in the interstitial is i^^lx = 20 Ry and for 
the potential EP°^^ = 196 Ry. For the (2 x 2) surface 
unit-cells, a (6 x 6 x 1) Monkhorst-Pack grid was used 
for the Brillouin zone integrations, for the gas-phase cal- 
culations in the rectangular supercells F-point sampling 



B. Stable adsorption sites and adsorption 
geometries 



In order to complement the preceding structure deter- 
mination by quantitative LEED Q a first series of cal- 
culations addressed the stable O and H adsorption sites 
by comparing the binding energy at the high-symmetry 
sites offered by the Ru(OOOl) surface. Consistent with 
previous calculations [l7| and in full agreement to the 
experimental analysis, the most favorable O adsorption 
site in purely oxygen-containing (2 x 2) overlayers with 
one up to four O atoms per surface unit-cell is the hep 
hollow site. On the contrary, for H atoms in the same 
(2 X 2) overlayers the most stable site is the fee hollow site, 
albeit with only a small and largely coverage- independent 
difference of about 50 meV/H atom to the hep site. Con- 
cerning the (2 X 2)-{0+E) and (2 x 2)-(0-F3H) coad- 
sorption phases, we subjected all structural models con- 
sidered in the recent quantitative LEED analysis of ref. 
to a full geometry optimization. The results are sum- 
marized in Table [J and are consistent with the preceding 
work in that the model with the lowest Pendry i?p-factor 
is also the one with the highest binding energy and the 
detailed geometric parameters of these two phases sum- 
marized in Table |TT] agree very well. This nicely confirms 
the assessment of ref. Q that LEED is indeed sensi- 
tive to the position of the H adsorbates. In these most 
stable structures the H atoms sit therefore, as already 
indicated in Fig. [U in the fee sites in the (2 x 2)-(0+II) 
configuration, and in the hep sites in the (2 x 2)-(0+3II) 
configuration. 
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TABLE I; Calculated total binding energy in meV and using gas-phase O2 and H2 as zero reference. Compared are the 
structural models for the (2 x 2)-(0+H) and (2 x 2)-(0+3H) coadsorption phases considered in ref. 0], using exactly the 
numbering employed in this preceding work. Model-X means that the structure relaxed to the structure of model-X. The most 
stable structure corresponding to the highest binding energy (marked in bold) is for both phases exactly the one identified in 
the quantitative LEED analysis of ref. [9|], and is shown in Fig. [l] 

1 2 3 4 5 6 7 8 

-Sb"ot"' 3400 3146 3301 model-2 model-3 model-2 model-1 3270 
-^b°ot™^ 3825 3049 3743 3752 2097 model-1 model-1 model-1 



TABLE IL Computed geometric parameters of the stable 
(2 X 2)-(0-|-H) and (2 x 2)-(0-|-3H) structures shown in Fig. [T] 
compared to the results of the experimental LEED structure 
determination of ref. [^. d is the adatom-Ru bond length, 
dij the mean vertical distance between layer i and j (with the 
adsorbate forming layer and the topmost Ru layer forming 
layer 1 etc.), and Zi is the vertical buckling of layer i. For 
the mean vertical layer distances, the values in brackets addi- 
tionally indicate the relative changes with respect to the bulk 
interlayer distance. 



Theory 



(2 x 2)-(0+H) 



Experiment 



d(H-Ru) (A) 1.89 1.83±0.13 
d(O-Ru) (A) 2.01 

Adoi(H-Ru) (A) 1.15 1.03-f0.15/-0.08 

Adoi(O-Ru) (A) 1.22 1.19±0.02 

Adi2 (A) 2.13 (-1.2%) 2.10±0.02 (-1.9 ±0.9%) 

Ad23 (A) 2.18 (+1.5%) 2.14±0.02 (0.0 ±0.9%) 

Ad34 (A) 2.17 (+1.0%) 

Azi (A) 0.1 0.1 

Az2 (A) 0.05 0.01 

Az3 (A) 0.03 



Theory 



(2 X 2)-(0+3H) 



Experiment 



d(H-Ru) (A) 1.86 - 1.88 1.89+0.3 
d(O-Ru) (A) 2.0 

Adoi(H-Ru) (A) 1.02 1.08+0.15 

Adoi(O-Ru) (A) 1.24 1.19+0.025 

Adi2 (A) 2.18 (+1.2%) 2.14+0.02(0.0 + 0.9%) 

Ad23 (A) 2.16 (+0.3%) 2.12+0.02 (-0.9 + 0.9%) 
Ad34 (A) 2.18 (+1.2%) 

Azi (A) 0.09 0.07 

Az2 (A) 0.02 0.03 
Az3 (A) 0.03 



IV. SCLS RESULTS AND DISCUSSION 

A. Comparison of computed and measured SCLSs 

From the geometric structures shown in Fig. [T] one ex- 
pects the foUowing first layer components in the mea- 
sured SCLS spectra: a peak due to Ru atoms without ad- 
sorbate coordination (Si) in the clean surface, as weh as 
peaks due to zero- and onefold 0-coordinated Ru atoms 
in the (2 x 2)-0 structure, Si and Si(0) respectively. 
In the (2 x 2)-(0+H) structure there are Ru atoms co- 



ordinated to 1 O + 1 Hfcc (Si(0+Hfcc)) and without 
adsorbate neighbor (Si); in the (2 x 2)-(0+3H) struc- 
ture there are Ru atoms coordinated to 1 O + 2 Hhcp 
(Si(0+2Hhcp)) and Ru atoms coordinated to 3 Hhcp 
(Si(3Hhcp); and in the (1 x 1)-H layer there are Ru atoms 
coordinated to 3 Hfcc neighbors (Si(3Htcc)). In addition 
there can be non-negligible second layer shifts (S2), where 
in the (2 x 2) layers one needs to furthermore distinguish 
between Ru atoms that are directly underneath an ad- 
sorbed O atom (S2(0)) and those that are not (S2). The 
SCLS data analysis and assignment for the clean surface 
and (2 x 2)-0 layer has already been established in ref. 
[3] and is the basis for the decomposition of the now mea- 
sured spectra shown in Fig. [3l Specifically, we followed 
the strategy explained in detail in ref. in which no 
attempt is made to resolve the small 82(0) second layer 
component in the 0-containing structure. Within this 
approach the assignment of the new peaks in the H over- 
layer and the two H-containing coadsorbate structures is 
straightforward. Also for the latter two only one non- 
negligible second layer peak was considered and we come 
back to the consequences for the derived SCLS values be- 
low. However, we note already here that although small, 
the shift of the first layer Si(0) component of the (2 x 2)- 
O overlayer when it becomes the Si(0+Hfcc) component 
in the (2 x 2)-(0+H) phase and even further when it is 
turned into the Si(0+2Hhcp) component in the (2 x 2)- 
(0+3H) phase is directly visible from the spectra and 
thus independent of the specific decomposition procedure 
employed. The same holds for the shift of the Si peak of 
the clean surface when turned into the Si(3Hfcc) compo- 
nent of the (1 X 1)-H layer, cf. Fig.d 

Table IIIII compiles the SCLS values derived from the 
decomposition shown in Fig. [3] and compares them to 
the calculated total shifts. The agreement is on the same 
remarkable level as that of ref. [J|, with maximum de- 
viations of the order of 40meV. As discussed in detail 
in this preceding work, one could presumably reach an 
even better agreement by explicitly considering the small 
82(0) component in the experimental spectra, which is, 
however, not necessary for the scope of this work. Essen- 
tially, the experiment-theory comparison is thus flawless, 
which reconfirms not only that the employed GGA-PBE 
functional seems to allow for a quantitative determina- 
tion of the quantity of interest to our study, but also fully 
supports the experimental peak assignment. On this ba- 
sis it is intriguing to realize that the derived first layer 
SCLS values are again quasi-linearly additive with re- 
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TABLE III: Comparison of measured and calculated SCLS values for the five investigated structures. All shifts are given 
in meV. See text for an explanation of the employed nomenclature. No attempt was made to fit a S2(0) component in the 
experimental spectra (see text). 
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TABLE IV: Calculated final-state contributions to the total 
SCLSs of first layer Ru atoms directly coordinated to zero, 
one, two or three adsorbed H atoms. 



1 2 

H coordination 



FIG. 4: Experimental H-induced changes of the first layer 
SCLSs as a function of the number of directly coordinated 
H atoms. Shown is the value by which the SCLS changes 
upon H neighbor addition, i.e. the reference is the Si value 
of the clean Ru(OOOl) surface for the purely H-coordinated 
Ru atoms, the Si value of the (2 x 2)-0 structure for the Ru 
atoms coordinated solely to H in the coadsorbate phases, and 
the Si(0) value of the (2 x 2)-0 structure for the Ru atoms 
coordinated to both O and H in the coadsorbate phases. 



spect to the coordination with H atoms. This is better 
visuahzed in Fig. 2] which shows the H-induced changes 
of the SCLSs, i.e. the amount by which a first layer Ru 
SCLS changes upon adding another H nearest neighbor. 
Strikingly, the observed additivity seems even indepen- 
dent of the presence of the O coadsorbate, and slightly 
different shifts are obtained for H neighbors in either fee 
or hep hollow sites. We furthermore mention that the 
O- and H-induced shifts are also closely the same to the 
ones observed in the (2 x 2)-(0-|-H) structure for a (2 x 2)- 
(0-l-H-l-CO) structure which additionally contains CO on 
top of those Ru surface atoms, which in the (2x2)-(0-|-H) 
have no direct adsorbate neighbors This suggests 

also a minor influence of near-next neighbor occupation 





Ofold 


ffold 2fold 3fold 


Clean 


-62 




(2 X 2)-0 


-45 




(2 X 2)-(0+Hfcc) 


-75 


-36 


(2 X 2)-(0-f3Hhcp) 




-38 -38 


(2 X 2)-4Htcc 




-44 



compared to the here discussed direct neighbor effects. 



B. Analysis of the "additivity rule" 

As a first step to analyze the physics governing the 
observed additivity rule it is useful to decompose the to- 
tal calculated shifts into initial- and final-state contri- 
butions. Table HVl compiles the final-state contributions 
to all first layer SCLSs of the five measured structures. 
Overall, the screening correction is found to be of the 
same order of magnitude as the actual total H-induced 
shifts presented in Fig. 21 as well as only marginally 
smaller than the screening contributions reported in ref. 
;4] for the 0-induced SCLSs at Ru(OOOl). This is dis- 
tinctly different to the situation for the initial-state shifts, 
which for the latter 0-containing structures were one 
order of magnitude larger than the now considered H- 
induced shifts. 

The consistently obtained negative sign of the screen- 
ing correction indicates that the screening capabilities of 
the core-ionized system are enhanced at the surface com- 
pared to the bulk. For the clean surface, the detailed 
analysis put forward in refs. [3, has already ratio- 
nalized this by the narrowing of the local valence c?-band 
induced by the lowered coordination of the surface atoms. 
This in turn leads in general to an enhancement of the 
density of states (DOS) at and above the Fermi level com- 
pared to the bulk situation, which — with the final-state 
effects at a transition metal surface like Ru(OOOl) largely 
due to the intra-atomic d-electron screening [19| — im- 
plies a more efficient screening at the surface. That this is 
not much affected by the presence of the adsorbates can 
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be understood considering the Anderson- Grimley-Newns 
model of chemisorption [20| : Consistent with this model 
the computed adsorbate induced DOS-change concerns 
primarily the formation of a bound state at the bot- 
tom of the substrate d-band and a broad (antibonding) 
resonance at the upper edge of the d-band. With the 
Fermi- level in Ru positioned at about 2/3 band- filling, 
this leaves the DOS at the Fermi-level and thereby the 
screening capabilities largely unaffected compared to the 
situation at the clean surface. In this understanding one 
would also not expect a strong influence of the adsor- 
bate induced geometric relaxation on the magnitude of 
the SCLS screening contribution. Corresponding calcu- 
lations in which the Ru atoms were frozen to their po- 
sitions in a bulk-truncated surface geometry and only 
the adsorbate positions were fully relaxed indeed confirm 
this expectation with computed final-state shifts that 
differ by less than lOmcV from the corresponding val- 
ues for the fully relaxed geometries summarized in Ta- 
ble IIVI Overall, the final-state contributions obtained 
for the 0-induced shifts and H-induced shifts are thus 
quite similar, with the values for different H coordina- 
tions showing even less variations. With a magnitude of 
~ — 55±20meV, consideration of these final-state effects 
is essential to reach the reported quantitative agreement 
with the measured total SCLSs. On the other hand, the 
scatter obtained in the final-state contributions for the 
different H coordination situations is too small to obscur 
the trend of about linearly increasing initial-state shifts 
with increasing H coordination, even though the latter 
are much smaller than the 0-induced ones where the ad- 
ditivity rule was observed before. 

As expected from the understanding developed in ref. 
Q we therefore arrive at the conclusion that the observed 
correlation of first layer SCLSs with the number of H 
nearest neighbors is an initial-state effect. The forma- 
tion of bonds between the adsorbate and the first layer 
Ru atoms affects the local valence d-band, the center of 
which shifts to maintain local charge neutrality. This 
changes the local (near nucleus) electrostatic field, which 
in turn acts on the core electrons as well and induces the 
finite SCLSs. Effectively, the about linear dependence 
of the latter on the number of H neighbors thus simply 
reveals that the type of bonding, i.e. the bonds formed 
per H nearest neighbor, is roughly the same in all over- 
layers studied. With all caveats given in ref. [3] this 
may be interpreted in the conceptual language of charge 
transfer in that the amount of charge transferred to each 
H atom remains approximately constant in the different 
structures. The much smaller electronegativity of H then 
also rationalizes the much smaller observed H-induced 
initial-state shifts (~ 50meV per H atom) compared to 
the 0-induced shifts reported before (~ 400 meV per O 
atom Q), and is fully consistent with the very different 
work function changes induced by the two adsorbates. In 
agreement to experiment (2T| we compute the latter to 
be essentially zero even for the full monolayer H cover- 
age, in stark contrast to measured |2^] and computed (l7j 



0-induced changes of up to ~ 1.5 eV for the full (1 x 1)- 
O monolayer. The thereby suggested near neutrality of 
H on Ru(OOOl) may also explain why the additivity rule 
prevails even in the denser adsorbate overlayers. For an 
effectively charged adsorbate, some degree of depolariza- 
tion, i.e. reduced charge transfer per adsorbate, would be 
favorable for higher coverages as it reduces the repulsive 
adsorbate-adsorbate interaction. This effect is indeed ob- 
served for dense O overlayers [1, [3| , but is too small there 
to significantly affect the near-linear SCLS correlation 
with O coordination. In this respect, one would expect 
a noticeable violation of the additivity rule particularly 
for adsorbates with only a modest bond strength (imply- 
ing smaller total magnitudes of the adsorbate-induced 
initial-state SCLSs) , but relatively strong induced dipole 
moment (suggesting more pronounced depolarization ef- 
fects at decreased adsorbate-adsorbate distances). This 
situation could for example be given for molecular ad- 
sorbates like CO or NO, and we are currently analyzing 
coadsorbate overlayers containing O, H and these molec- 
ular species in this regard p^ . 



With the SCLSs sensitively reflecting the formed 
adsorbate-substrate bonds, it is clear that they can also 
be used as local probes of the adsorption site. This was 
already shown in the preceding work on the 0-containing 
structures, where the computed Ru SCLSs for structures 
with O atoms in fee sites were consistently about 100-200 
meV smaller than the corresponding values in structures 
with O atoms in hep sites and therewith irreconcilable 
with the experimental data For O adsorption into 
hep sites the induced SCLS change per O neighbor is thus 
slightly larger, whereas, interestingly, for Rh(lll) where 
the fee site is the more stable adsorption site the situa- 
tion is exactly reversed in that there O adsorption into fee 
sites leads to slightly larger SCLS increases [3] . For both 
surfaces, a somewhat more "negatively charged" oxygen 
atom is thus adsorbed in the more stable adsorption site, 
suggesting a preference for a stronger ionic bonding. In 
this respect, it is intriguing to notice from Fig.|4]that also 
for the H overlayers the slightly more stable fee site is the 
one that leads to slightly larger induced SCLS changes. 
The less favorable hep site is only populated in the O- 
containing dense (2 x 2)-(0-|-3H) coadsorbate structure 
and it is tempting to attribute this stabilization in parts 
to that a slightly less ionic hydrogen is preferred at the 
corresponding short inter-adsorbate distances. While in 
this structure the effect concerns only the position of the 
H atoms with the O atoms remaining in the hep sites, 
coadsorption of CO is able to shift also the O atoms to 
fee sites. Analyzing the correlation of the respective sta- 
bilities of different sites with the induced SCLS changes 
is thus another aspect of the mentioned on-going work 
on coadsorption structures containing O, H, CO and NO 

0. 
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V. CONCLUSIONS 

In conclusion we have presented a joint experimental 
and theoretical study of Ru(OOOl) surface core level shifts 
induced by O and H adsorbates. Remarkable quanti- 
tative agreement concerning the energetics, geometries, 
and SCLS values of the investigated (1 x 1)-H, (2 x 2)- 
(O+H), and (2 x 2)-(0+3H) (co)adsorbate overlayers on 
Ru(OOOl) is achieved. The quasi-linear dependence of 
the adsorbate-induced SCLSs on the number of directly 
coordinated nearest neighbors reported before for O over- 
layers on Ru(OOOl) H, Rh(lll) [aj and Rh(lOO) 0, as 
well as for H adsorbate layers on all Rh low-index sin- 
gle crystal surfaces [5j] is found to prevail also in the H- 
(co)adsorbate structures on Ru(OOOl). The computed 
decomposition of the total SCLSs into initial- and final- 
state contributions reveals an enhanced screening at the 
surface of the (co) adsorbate layers that is of the same 
order of magnitude as for the pure O overlayers studied 
before. Accounting for this final-state correction is thus 
necessary to obtain the reported quantitative agreement 
with the measured data, while its actual value exhibits 
only a small scatter for the differently coordinated surface 
atoms. In accordance with the detailed analysis put for- 
ward in ref. [3l , the observed additivity of the adsorbate- 
induced SCLSs is thus a pure initial-state effect. In terms 
of the conceptual charge transfer language it merely re- 
flects a roughly contant amount of charge transferred to 
each directly coordinated adsorbate neighbor. Not sur- 
prisingly, this amount of charge (and concomitantly the 
magnitude of the observed SCLS change) is significantly 
smaller for the almost electroneutral H adsorbates than 
for the electronegative O adsorbates. Intriguingly, the 
SCLS data reveals slightly larger induced changes for H 
atoms adsorbed in fee sites than in hep sites, consistent 
with the interpretation that a slightly more ionic H ad- 
sorbs in the prior sites. The established sensitivity even 
to such highly similar adsorption sites, as well as the 
prevalence of the additivity rule confirms the suggested 
use of SCLSs as directly accessible and valuable finger- 
prints of the adsorbate populations. 

Appendix 

Similar to the analysis performed in ref. [3| , the exper- 
imental data were fitted with a combination of Doniach- 
Simjic functions [11] convoluted with a Gaussian broad- 
ening and a linear background. The lineshape parame- 



ters are the Lorentzian (L), Gaussian (G) and asymmetry 
(a). The following outhnes the fitting procedure used to 
extract these parameters and the values of the Ru 8^5/2 
SCLSs. 

The absolute binding energy of the peaks was cali- 
brated to the Fermi-level, measured for each spectrum. 
Moreover, the intensity of the spectra was normalized to 
the background at the low binding energy side (279.5 eV) 
in order to account for variations of the photon beam 
intensity. The best fit parameters were extracted from 
the fit of all the data taken at the three photon energies 
used. First the spectrum of the clean surface was fitted 
with three components: bulk {L = 0.20 eV, G = 0.06 eV, 
a = 0.07), Si (L = 0.24eV, G = 0.06eV, a = 0.05), 
and S2 {L = 0.20 cV, G = 0.06 eV, a = 0.08). The ab- 
solute position of the bulk peak was determined as well 
and this position was kept fixed in all other fits. Next, 
we fitted the (2 x 2)-0 spectrum to determine the Si(0) 
lineshape (L = 0.24 eV, G = 0.12 eV, a = 0.06), which 
is wider than the Si component of the clean surface as 
already found in our previous work The Si peak in- 
tensity turns out to be ~ 1/4 of the same peak of the 
clean surface, as it should be at the oxygen coverage of 
0.25 monolayer. Moreover, the intensity of the bulk peak 
is increased while that of S2 is lowered, consistent with 
the expectation that a near-zero and unresolved S2(0) 
peak contributes to the bulk region. 

We then continued with the (1 x 1)-H structure where 
again three components are expected: bulk, S2 (same 
lineshape parameters as S2 in the clean surface) and 
Si(3Hfcc) {L = 0.24eV, G = O.lOeV, a = 0.05). It 
has to be noted here that the intensity of the different 
components is almost the same as for the clean surface, 
because the hydrogen atoms present on the surface have 
only small scattering effects on the electrons photoemit- 
ted from the Ru substrate. On the other hand, the pres- 
ence of H at the surface seems to affect the peak shape 
in the spectra, and we obtain significantly improved fits 
when allowing for a larger Gaussian width. As mentioned 
in the text, we interpret this by inhomogeneous broaden- 
ing due to a less than perfect ordering achieved for the H 
overlayer. The same situation (unchanged intensity with 
respect to the (2 x 2)-0, increased Gaussian width) is 
found for the coadsorbate structures, and the line shape 
parameters determined are: Si(O-l-lHfcc) {L = 0.24 eV, 
G = 0.15 eV, a = 0.06) in the (2 x 2)-(0-HH) structure; 
and Si(0+2Hhcp) [L = 0.26 eV, G = 0.16 eV, a = 0.06) 
and Si(3Hhcp) [L = 0.24 eV, G = 0.12 eV, a = 0.05) in 
the (2 X 2)-(0+3H) structure. 
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